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ABSTRACT 

A four-dimensional statistical description of electromagnetic radiation is developed and applied to the 
analysis of radio pulsar polarization. The new formalism provides an elementary statistical explanation 
of the modal broadening phenomenon in single pulse observations. It is also used to argue that the degree 
of polarization of giant pulses has been poorly defined in past studies. Single and giant pulse polarimetry 
typically involves sources with large flux densities and observations with high time resolution, factors 
that necessitate consideration of source-intrinsic noise and small-number statistics. Self noise is shown 
to fully explain the excess polarization dispersion previously noted in single pulse observations of bright 
pulsars, obviating the need for additional randomly polarized radiation. Rather, these observations are 
more simply interpreted as an incoherent sum of covariant, orthogonal, partially polarized modes. Based 
on this premise, the four-dimensional covariance matrix of the Stokes parameters may be used to derive 
mode-separated pulse profiles without any assumptions about the intrinsic degrees of mode polarization. 
Finally, utilizing the small-number statistics of the Stokes parameters, it is established that the degree 
of polarization of an unresolved pulse is fundamentally undefined; therefore, previous claims of highly 
polarized giant pulses are unsubstantiated. 



Subject headings: methods: data analysis 
techniques: polarimetric 

1. INTRODUCTION 



methods: statistical — polarization — pulsars: general 



Radio pulsars exhibit dramatic fluctuations in total and 
polarized flux densities on a diverse range of longitudi- 
nal, temporal and spectral scales. As a function of pulse 
phase, both average profiles and sub-pulses make sud- 
den transitions between orthogonally pola rized modes (e.g. 
iTavlor et aLlll971l : [Manchester et al.lll975ft . The radiation 
at a single pulse phase often appears as an incoherent su- 
p erp osition of modes, both orthogonal and non- orthogonal 
fe.g. lBacker fc Rankidfl980l lMcKinnonll2003aD . Evidence 
has also been presented for variations that may indicate 
stochastic g eneralized Faraday rotation in the pulsar mag- 
netosphere (|Edwards fc Stappersll2004) . 

In an extensive stud y of single-pulse polariz ation fluc- 
tuations at 1400 MHz, IStinebring et all (|1984D observed 
that histograms of the linear polarization position angle 
were broader than could be explained by instrumental 
noise alone. This modal broadening of the position angle 
distribution was interpreted as e vidence of a superposed 
rando mized emission component. iMcKinnon fc Stinebrina 
(1998) revisited the issue with a statistical model that de- 
scribed correlated intensity fluctuations of completely po- 
larized orthogonal modes. Although the predicted position 
angle distributions were qualitatively similar to the obser- 
vations, the measured histograms were wider than those 
produced by a numerical simulation of modal broadening 
that included instrumental noise. Source-intrinsic noise 
was later used to explain the excess polarization scatter 
of bright pulses in an analysis of mode-separated profiles 
(iMcK innon fc Stinebringll2000t ): however, it was given no 
further c onsideration in a l l subsequent statis tical treat- 
ments bv lMcKinnoil (12001 l2003bl I2004L I2006T) . 

These begin with the reasonable assumptions that the 



instantaneous signal-to-noise ratio S/N is low and a suffi- 
ciently large number of samples have been averaged, such 
that the stochastic noise in all four Stokes parameters 
can be treated as uncorrelated and normally distributed. 
These assumptions form part of the three-dimensional 
eigenvalue analys es of the Stokes polarization vect or by 



Mc Kinnonl (l2004h a nd lEdwards fc Stappersl (|2004h . As 



in IStinebring et al] (|1984h , these studies concluded that 
modal broadening is due to the incoherent addition of 
randomly polarized radiation that is intrinsic to the pul- 
sar. The basic premises of these experiments are valid for 
the vast majority of pulsar observations; however, they 
become untenable for the bright sources on which these 
studies focus. That is, when the instantaneous S/N > 1, 
the Stokes polarization vector can no longer be treated in 
isolation and the correlated self noise intrinsic to all four 
stochastic Stokes parameters must be accounted. 

These considerations are particularly relevant to the 
study of giant pulses. Those from the Crab pulsar reach 
brightness temperatures in excess of 10 37 Kelvin and re- 
main urn^sol_yed_j^ol2SCTvations with nanosecond resolu- 
tion (jHankins et al.l 120 03). Consequently, previous anal- 
yses of giant pulses have typically presented polarization 
data at the sampling resolution of the instrument; that is, 
where the ti me-bandwidth prod uct is of the order of unity. 
For example, Hcil es et al.l (|1970f ) studied giant pulses from 
the Crab pulsar using an instrument with bandwidth, 
8v = 8.3 kHz, and presented plots of the Stokes pa- 
ramet ers at a resolution of t = 120^s. ICognard et alj 
( 199(f) plotted the total intensity and circular polarization 
of the strongest giant pulse and interpulse observed from 
PSR B1937+21 with r = 1.2^s and 5v = 5 00 kHz. Using 
a base band recorder with 5v = 500 MHz, lHankins et al.l 
(2003|) plotted the intensities of left and right circularly 
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polarized radiation from Crab giant pulses at a time reso- 
lution of 2 ns. Each of these studies concluded that giant 
pulses are highly polarized. However, in each experiment, 
t5v *~ 1; at this resolution, every discrete sample of the 
electric field is completely polarized, regardless of the in- 
trinsic degree of polarization of the source. That is, the 
instantaneous degree of polarization is fundamentally un- 
defined. 

To study the polarization of intense sources of impulsive 
radiation at high time resolution, it is necessary to con- 
sider averages over small numbers of samples and source- 
intrinsic noise statistics. These limitations are given care- 
ful at tention in the statistical theory of polarized shot 
noise ([Corded 119761) . This seminal work enables charac- 
terization of the timescales of microstructure polarization 
fluctuations via the auto-correlation functions of the total 
inte nsity and degrees of line ar and circular polarization 
(e.g. ICordes fc Hankind [19771 ). It has also been extended 
to study the degree of polarization via the cross-correlation 
as a function of time lag betwee n the instantaneous total 
intensity spectra of giant pulses (|Cordes et al.l [20041 ). 

This paper presents a complimentary approach based on 
the four-dimensional joint distribution and covariance ma- 
trix of the Stokes parameters, with emphasis placed on the 
statistical degrees of freedom of the underlying stochas- 
tic process. Relevant theoretical results are drawn from 
various sources, ranging from studies of the scattering of 
mono chromatic light in optical fibres (e.g. ISteeger et alj 
1 19841 ) to the classification of sy nthetic aperture radar im- 
ages (e.g. iTouzi fc Lopes1ll996l ). 

Following a brief review of polarization algebra in § [2] 
the joint probability density functions of the Stokes param- 
eters at different resolutions are derived in § [31 where the 
results are compared and contrasted with previous works. 
The formalism is related to the study of radio pulsar polar- 
ization in § [4[ where the effects of amplitude modulation 
and wave coherence on the degrees of freedom of the pul- 
sar signal are discussed. Finally, the results are utilized to 
reexamine past statistical analyses of orthogonally polar- 
ized modes, randomly polarized radiation, and giant pulse 
polarimetry. It is concluded that randomly polarized ra- 
diation is unnecessary and the degree of polarization of 
giant pulses must be more rigorously defined. Potential 
applications of the four-dimensional statistics of polarized 
noise are proposed in § [5) 

2. POLARIZATION ALGEBRA 

This section reviews the relevant algebra of both Jones 
and Mueller representations of polarimetric transforma- 
tions. Unless otherwise noted, the notation and terminol- 
ogy synthesizes that of th e similar appro a ches to polariza - 
tion algebra present ed by ICloudd (|1986l ). iBrittonl (|2000D . 
and lHamakel (|2000D . 

2.1. Jones Calculus 

The polarization of electromagnetic radiation is de- 
scribed by the second-order statistics of the transverse 
electric field vector, e, as repre sented by the compl ex 2x2 
coherency matrix, p = (ee') (|Born fc Wolfl [T980). Here, 
the angular brackets denote an ensemble average, is 
the Hermitian transpose of e, and an outer product is 
implied by the treatment of e as a column vector. A 



useful geometric relationship between the complex two- 
dimensional space of the coherency matrix and the real 
four-dimensional space of the Stokes parameters is ex- 
pressed by the following pair of equations: 

p = S k a k /2 (1) 
S fe = Tr(<7 fe p). (2) 
Here, S k are the four Stokes parameters, Einstein notation 
is used to imply a sum over repeated indeces, < k < 3, 
(Tq is the 2x2 identity matrix, cri-3 are the Pauli ma- 
trices, and Tr is the matrix trace operator. The Stokes 
four-vector is composed of the the total intensity So and 
the polarization vector, S — (Si, S2, S3). Equation (fT]) 
expresses the coherency matrix as a linear combination 
of Hermitian basis matrices; equation ([2]) represents the 
Stokes parameters as the projections of the coherency ma- 
trix onto the basis matrices. These properties are used 
to interpret the well-studied statistics of random matrices 
through the familiar Stokes parameters. 

Linear transformations of the electric field vector are 
represented using complex 2x2 Jones matrices. Substitu- 
tion of e! — Je into the definition of the coherency matrix 
yields the congruence transformation, 

p' = Jpjt, (3) 
which forms the basis of the various coordinate transfor- 
mations that are exploited throughout this work. If J is 
non-singular, it can be decomposed into the product of a 
Hermitian matrix and a unitary matrix known as its polar 
decomposition, 

J = JB*(i8)Rft(^ (4) 

where J = IJI 1 ^ 2 , B T f 1 ( / 3) is positive-definite Hermitian, 
and Rfj(0) is unitary; both Bm(/5) and R^(0) are uni- 
modular. Under the congruence transformation of the co- 
herency matrix, the Hermitian matrix 

B r ?,, ((3) — exp(/3 rh • cr) = cr cosh (3 + rh • cr sinh (3 (5) 
effects a Lorentz boost of the Stokes four- vector along the 
rh axis by a hyperbolic angle 2(3. As the Lorentz transfor- 
mation of a spacetime event mixes temporal and spatial 
dimensions, the polarimetric boost mixes total and polar- 
ized intensities, thereby altering the degree of polarization. 
In contrast, the unitary matrix 

Rfi {4>) = exp(z0 n ■ cr) — er cos (f> + in • cr sin 4> (6) 
rotates the Stokes polarization vector about the n axis by 
an angle 2cj). As the orthogonal transformation of a vector 
in Euclidean space preserves its length, the polarimetric 
rotation leaves the degree of polarization unchanged. 

These geometric interpretations promote a more in- 
tuitive treatment of the matrix equations that typically 
arise in polarimetry Boost transformations can be uti- 
lized to convert unpolarized radiation into partially polar- 
ized radiation, and rotation transformations can be used 
to choose the orthonormal basis that maximizes symme- 
try. These properties are exploited in § [3] to simplify the 
relevant mathematical expressions that describe the four- 
dimensional joint distribution of the Stokes parameters. 

2.2. Eigen Decomposition 

It proves useful in § [3] to express the coherency matrix 
as a similarity transformation known as its eigen decom- 
position, 

' = AO*" ^ 
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Here, R = (eo ei) is a 2 x 2 matrix with columns equal to 
the eigenvectors of p, and A m are the corresponding eigen- 
values, given by A = (5 ± |S|)/2 = (1 ± P)S /2, where 
P = Sq/\S\ is the degree of polarization. If the signal is 
completely polarized, then \i = 0. If the signal is unpo- 
larized, then there is a single 2-fold degenerate eigenvalue, 
A = So/ 2 and R is undefined. 

If the eigenvectors are normalized such that e k 'e k = 1, 
then equation ([7]) is equivalent to a congruence transforma- 
tion by the unitary matrix, R. In the natural basis defined 
by IV , the eigenvalues A m are equal to the variances of two 
uncorrelated signals received by orthogonally polarized re- 
ceptors described by the eigenvectors. The total intensity, 
So = Ao + Ai; the polarized intensity, S± — \S\ = Ao — Ai; 
and 52 = S3 — 0. That is, Rt rotates the basis such that 
the mean polarization vector points along Si, providing 
cylindrical symmetry about this axis. 

2.3. Mueller Calculus 

The congruence transformation of the coherency matrix 
by any Jones matrix J may be represented by an equiv- 
alent linear transformation of the Stokes parameters by a 
real-valued 4x4 Mueller matrix 



1 



M^-Tr^Jo-fcjt), 



such that 



S : 



M*S k . 



(8) 



(9) 



Mueller matrices that have an equivalent Jones matrix are 
called pure, and suc h transformations are related to the 
Lorentz group (e.g. [Barakat 19631: iBrittoi] I2000D . This 
motivates the definition of an inner product 



(A, B) = A k B k = V kk A k B k = A B A • B, (10) 

where A and B are Stokes four-vectors and rjij is the 
Minkowski metric tensor with signature (+,—,—,—). The 
Lorentz invariant of a Stokes four-vector is equal to four 
times the determinant of the coherency matrix; that is, 

5 2 ^<5,5) = 5 2 -|S| 2 = 4|p| (11) 
Similarly, the Euclidean norm ||5|| is twice the Frobenius 
norm of the coherency matrix ||p||; i.e. 

||5|| 2 ^5 2 + |S| 2 =4||p|| 2 . (12) 

The coherency matrix is a positive semi-definite Hermitian 
matrix; therefore, the Lorentz invariant of any physically 
realizable source of radiation is greater than or equal to 
zero. It is equal to zero only for completely polarized radi- 
ation, when the degree of polarization is unity (So — \S\). 
As with the spacetime null interval, no linear transforma- 
tion of the electric field can alter the degree of polarization 
of a completely polarized source. 

3. JOINT DISTRIBUTION FUNCTION 

In this section, the joint distribution functions of the 
Stokes parameters for a stationary stochastic source of po- 
larized radiation are derived. Three regimes of interest 
are considered: single samples, local means, and ensemble 
averages of large numbers of samples. The electric field 
vector is assumed to have a jointly n ormal density func- 
tion as described bv I Goodman! (|1963[ ). This distribution 
may not accurately describe the possibly non-linear elec- 
tric field of intense non-thermal radiation. For example, 



over a limited range of pulse phase, the fluctuations in the 
intensity of the Vela pulsar have a lognormal distribution 
that is consistent with the predictions of stochastic growth 
theory (|Cairns et alj 12001). The normal distribution is a 
valid choice if the pulsar signal can be accurately mod- 
eled as polarized shot-n oise, provided that the density of 
shots is sufficiently high (|Cordeslll976l ). It is also the min- 
imal assumption if no information about the higher-order 
moments of the field is available. 

3.1. Single Samples 

As second-order moments of the electric field, the co- 
herency matrix and Stokes parameters are defined only 
after an average is made over some number of samples. 
Given a single instance of the electric field e (for example, 
discretely sampled at an infinitesimal moment in time) 
define the instantaneous coherency matrix, r = ee\ and 
Stokes parameters, s k — Tr(cr fe r), such that the ensemble 
averages, p = (r) and S k = (s k ). It is trivial to show 
that the determinant of the instantaneous coherency ma- 
trix, |r| = 0, regardless of the degree of polarization of the 
source or the probability density of the electric field. That 
is, the instantaneous degree of polarization is undefined. 

The complex-valued components of e are the analytic 
signals associated with the real-valued voltages measured 
in each receptor. If the voltages are normally distributed 
with zero m ean, then e has a bivariate complex normal 
distribution (| Goodman! 1 19631 ). 

/(e) = ^exp(-etp- 1 e), (13) 



ir 2 \p\ 



where p = S k <r k /2 is the population mean coherency ma- 
trix. To compute the probability density function of the 
instantaneous Stokes parameters, note that /(e) is inde- 
pendent of the absolute phase of e. Conversion to polar 
coordinates and marginalization over this variable yields 
the intermediate result, 

/(|e |, | ei U) = exp (- e tp-i e ) , (14) 



TT\P\ 



where tp is the instantaneous phase difference between eo 
a nd e\. Th i s resu lt may be compared with equation (2.10) 
of iBarakatl (|1987l ). except that here the evaluation of the 
inner product has been postponed. The three remaining 
degrees of freedom are described by the instantaneous po- 
larization vector s, for which the Jacobian determinant 



J(s; |e |, \ei\,ip) 



ds 



Application of the above with 
P 1 = TTTTT^o c - S 



8|eo|| ei || S |. (15) 



2P 



<t) = ^S k a k 



and 



e ] p~ x e = Tr (e f p~ 



l e =Tr p" 



(16) 



(17) 



yields the four-dimensional joint distribution of the instan- 
taneous Stokes parameters, 



f(s) 



S(s -\s\) 



exp 



2(5, a) 



(18) 



nS^so * V S 2 
where S is the Dirac delta function and so is the instanta- 
neous total intensity . This result is consistent with equa- 
tion (11) of lElivah u (1994). Converting to spherical polar 
coordinates, the inner product (5, s) = so(So — |S|)cos0, 
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where 6 is the angle between s and S. Subsequent in- 
tegration over s yields the marginal distribution of the 
instantaneous total intensity, 

2 / 2s S \ . , f 2s \S\ \ 
f( s o) ~ 757 ex P ( 55- smh —55— ' ( 1Q ) 



|s| ~*" V s 2 J V s 2 J ' 

which has mean So and variance 115*11 2 /2. Th i s dist ribu- 
tion is consisten t with equation (13 ) oflMandell (|l963l ) and 
equation (4) of IFercher fc Steegerl (|l98lh . As noted by 
these authors, the total intensity becomes % 2 distributed 
in the two cases of unpolarized and completely polarized 
radiation. This is easily seen in the natural basis defined 
by the eigen decomposition of the coherency matrix. As 
the squared norm of a complex number, the instantaneous 
intensity in each of the two orthogonally polarized modes 
is x 2 distributed with two degrees of freedom. If the sig- 
nal is unpolarized, then each mode contributes identically 
and the total intensity is \ 2 distributed with four degrees 
of freedom. If the signal is 100% polarized, then only one 
mode contributes and the total intensity is \ 2 distributed 
with two degrees of freedom. 

The marginal distributions of the instantaneous Stokes 
polarization vector components are most easily derived in 
the natural basis, where {S, s) = sqSq — siS±. Convert- 
ing s to cylindrical coordinates with the axis of symmetry 
along Si then integrating equation ([IS]) over the radial and 
azimuthal dimensions (as well as the total intensity) yields 
the marginal distribution of the instantaneous major po- 
larization, 



/Oi) = Trexp 



2\ Sl \S \ {2 Sl \S\ 
exp 



2 . . 2 f . (20) 

This is an asymmetric Laplace density with mean S\ 
and variance ||5| | 2 /2, c onsistent with equation (8) of 
IFercher fc Steegerl (|l98lh . The marginal distribution of 
the instantaneous minor polarization S2 is derived in Ap- 
pendix [5] by symmetry, equation (|A4|) yields 



/(S2,a) = 4 exp 



2|«2, 3 | 



(21) 



where S is the Lorentz interval. This is a symmetric 
Laplace density with mean zero and vari ance S 2 /2, as in 
equation (16) of IFercher fc Steegerl (|198lD . 

Equations (TiU[) through (f2"Tj) are plotted in Figure [1] 
The bottom two panels illustrate the asymmetric, three- 
dimensional Laplace distribution of the instantaneous po- 
larization vector s; the top panel shows the distribution 
of the magnitude of this vector sq = \s\. Qualitatively, in- 
stances of s are distributed as a tapered needle that points 
along the major polarization axis, with the greatest den- 
sity of instances in the head of the needle at the origin. 
For unpolarized radiation, the distribution is spherically 
symmetric and /(si) = /(s2,3)- For completely polarized 
radiation, /(s2,3) = ^(^2,3) and s± — sq; that is, the dis- 
tribution of s becomes one-dimensional and exponentially 
distributed along the positive Si axis. 

3.2. Local Means 

The distribution of the local mean of n > 1 independent 
and identically distributed instances of the Stokes param- 
eters s are derived from the distribution of the local mean 
coherency matrix 



1 



e,;ej 



1 



(22) 



which has a complex Wishart distribution with n degrees 
of fre edom (|Wishardll92l^lGoodmanlll963l ; lTouzi fc Loped 
1996); i.e. 



m = 



n 2n \p\ n ~ 2 



■exp(-nTr . (23) 



7TT(n)T(n~ l)\p\' 

The Wishart distribution is a multivariate generalization 
of the x 2 distribution, and plays an important role in com- 
munications engineering and information theory. As in 
§ 13.11 the Jacobian determinant 



J(S k ,p) = 



dSu 



dp 



= 8 



(24) 



and equation (|16[) are used to arrive at the joint distribu 
tion function of the local mean Stokes parameters, 

2n 2n^2(n-2) ( 2n(S,S 



7rr(n)r(n- l)^ 2 ™' 6 ^ V S 2 
In Appendix [B] this joint distribution is used to derive 
the probability density and the first two moments of the 
local mean degree of polarization, p — \s\/so- In Figure [21 
the distribution of p is plotted as a function of the intrin- 
sic degree of polarization P and the number of degrees of 
freedom n. Note that, when P = 1, f(p) = 5(p— 1); there- 
fore, this case is not shown. The expected value of p as 
a function P and n is plotted in Figure [3l the standard 
deviation o~ p is similarly plotted in Figure 3J 

Figure [3] demonstrates that the self noise intrinsic to a 
stationary stochastic source of polarized radiation induces 
a bias, (p) — P, in the estimated degree of polarization. The 
bias and standard deviation define the minimum sample 
size required to estimate the degree of polarization to a 
certain level of accuracy and precision. Given the sample 
size and a measurement of p, the probability density of p 
(eqs. |B2j and |B3| ) or the expectation value of p (eq. |B4j ) 
could be used to numerically estimate the intrinsic degree 
of polarization using metho ds similar to those reviewed by 
ISimmons fc Stewarti (|1985T ). 

3.3. Ensemble Averages 

By the central limit theorem, at large n the mean Stokes 
parameters tend toward a multivariate normal distribution 

1 ( 1 



(25) 



f(S) 



47T 2 |C| 1 /2 



exp 



-AS T C- 1 AS 



(26) 



where AS = S — S and C is the covariance matrix of 
the Stokes parameters. To derive the components of C, 
consider a completely unpolarized, dimcnsionless signal 
with unit intensity (that is, with mean coherency matrix 
Pu = CT o/2 and mean Stokes parameters S u = [1,0,0,0]) 
and covariance matrix C u = ^ 2 I, where ^ 2 is the dimen- 
sionless variance of each Stokes parameter and I is the 
4x4 identity matrix. This unpolarized signal may be 
transformed into any partially polarized signal with mean 
coherency matrix p via a congruence transformation 

p = bp^ = bbt 12 = b 2 /2, (27) 

wher e b is the Hermitian square root of 2p (e.g. iHamakerl 
2000) and the elements of p have physical dimensions 
such as flux density. The covariance matrix of the Stokes 
parameters of the resulting partially polarized signal is 
C = BC U B T , where B is the Mueller matrix of b, defined 
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Instantaneous Intensity, s /S 




-1.5 -1 -0.5 0.5 1 1.5 



Instantaneous Minor Polarization, s 33 /S 

Fig. 1. — Probability densities of the instantaneous Stokes parameters as a function of the degree of polarization, P. From top to bottom 
are the density functions of the total intensity, the major polarization, and the minor polarizations. In each panel, the unpolarized case is 
labeled P = 0, followed by P = 0.5, P = 0.8 and, where applicable, P = 1.0. 
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Fig. 2. — Probability densities of the local degree of polarization as function of the intrinsic degree of polarization, P, and the number of 
degrees of freedom, n. In each panel, the unpolarized case is labeled P = 0, followed by P = 0.5 and P = 0.8. 
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Fig. 4. — Standard deviation of the local degree of polarization as a function of the intrinsic degree of polarization. The number of degrees 
of freedom is printed above each curve. 
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by equation ([8]). Noting that B is symmetric, C = <^ 2 B 2 
is simply a scalar multiple of the Mueller matrix of p; i.e. 

dj = <r 2 (2SiSj - rjijS 2 ) . (28) 
This result is a genera lization of equation (9) in 
IBrosseau fc Barakatj ( 1992fl . which derives from the com- 
plex Gaussian moment theorem. In contrast, equation (|28[) 
requires no assumptions about the distribution of the elec- 
tric field. 

The covariance matrix of the Stokes parameters has its 
simplest form in the natural basis defined by the eigen 
decomposition of the coherency matrix, where 



a = e 




(29) 



Comparison between the diagonal of C and the variances 
of the distributions derived in § 13.11 yields the relation- 
ship between the dimcnsionlcss variance and the number 
of degrees of freedom, s 2 = (2n) _1 . In the natural basis, 
it is also readily observed that the polarization vector S 
is normally distributed as a prolate spheroid with axial 
ratio, e = ((1 + P 2 )/(l - P 2 )) 1/2 . The dimension of the 
major axis of the spheroid is equal to the standard devia- 
tion of the total intensity and the dimension of the minor 
axis goes to zero as the degree of polarization approaches 
unity. Furthermore, the multiple correlation between So 
and S is 

2P 

R = TT^- (30) 

The multiple correlation ranges from for an unpolarized 
source to 1 for a completely polarized source. It character- 
izes the correlation between total and polarized intensities 
and expresses the fact that the stochastic Stokes parame- 
ters cannot be treated in isolation. 

4. APPLICATION TO RADIO PULSAR ASTRONOMY 

The previous section develops the four-dimensional 
statistics of the Stokes parameters intrinsic to a single, 
stationary source of stochastic electromagnetic radiation. 
To apply these results to radio pulsar polarimetry, it is 
necessary to consider various observed properties of pul- 
sar signals, including amplitude modulation, wave co- 
herence, and the superposition of signals from multiple 
sources, such as instrumental noise and orthogonally polar- 
ized modes. Emphasis is placed on the statistical degrees 
of freedom of the radiation; in particular, the effects of am- 
plitude modulation and wave coherence on the identically 
distributed and independent conditions of the central limit 
theorem are considered. Attention is then focused on two 
areas of concern: randomly polarized radiation and giant 
pulse polarimetry. 

4.1. Amplitude Modulation 

Radio pulsar si gnals are wel l described as amplitude- 
modulated noise (jRickettj Il975f) . The modulation index 
is generally defined as (3 = ay/ So, where do and So are 
the standard deviation and mean of the total intensity. In 
single-pulse observations of radio pulsars, <7o is typically 
corrected for the instrumental noise estimated from the 
off-pulse baseline and the self noise is assumed to be negli- 
gible. A value of f3 greater than zero is then interpreted as 



evidence of amplitude modulation. However, for intense 
sources of radiation, the self noise of the ensemble average 
total intensity, a = s\\S\\ (cf. eq. [29] of Cordes 1976) 
must be taken into account. That is, the self noise of the 
source induces a positive bias in the modulation index. 

Amplitude modulation modifies the statistics of all four 
Stokes parameters. For example, under scalar amplitude 
modulation, the covariance matrix of the Stokes parame- 
ters becomes (A 2 )C, where A is the instantaneous inten- 
sity of the modulating function and (A 2 ) > (A) 2 = 1. 
That is, amplitude modulation uniformly increases the co- 
variances of the Stokes parameters. This can be inter- 
preted as a reduction in the statistical degrees of freedom 
of the signal, which becomes dominated by the fraction of 
realizations that occur when A ^> 1. That is, especially 
when the modulations are deep, the samples are not iden- 
tically distributed and the central limit theorem does not 
trivially apply. 

4.2. Wave Coherence 

The degrees of freedom of a stochastic process are also 
reduced when the samples are not independent. For plane- 
propagating electromagnetic radiation, statistical depen- 
dence is manifest in various forms of wave coherence, in- 
cluding that between the orthogonal components of the 
wave vector (i.e. polarization) and that between instances 
of the field at different coordinates (e.g. spectral, spa- 
tial, temporal). In radio pulsar observations, wave co- 
herence properties may be modified by propagation in 
the pulsar magnetosphe re and the interstellar medium 
(e.g.lLvutikov fc Parikbll2000t iM drose fc Macquart] |1998l; 
iMacauart fe Melrosdl200Ci) . 

Given the coherence time r c (|Mandell 119581 11959( ) of a 
band-limited source of Gaussian noise, the effective num- 
ber of degrees of freedom is 

T T 

N cS = — <-=n, (31) 

. T c + T T 

where T is the integration interval, r is the sampling in- 
terval, and n is the number of discrete time samples in 
the interval T. That is, owing to wave coherence, the sig- 
nal could be encoded by N e f{ independent samples without 
any loss of information. 

Substituting n = N e g into equation (|25[) or <; 2 oc 1/N e ff 



into equation (|26[) . it is seen that wave coherence inflates 
the four-dimensional volume occupied by the joint distri- 
bution of the mean Stokes parameters. With respect to the 
statistics of independent samples, this inflation increases 
both the modulation index of the mean total intensity and 
the eigenvalues of t he covariance matrix of the mean Stokes 
polar ization vector (|Edwards fc Stap pcrs 2004; McKi nnonl 
2004). Referring to Figure [31 it is readily observed that 
wave coherence also increases the local mean degree of po- 
larization. 

If only the covariance matrix of the Stokes parame- 
ters is measured, the effects of wave coherence are in- 
distinguishable from those of amplitude modulation. In 
fact, the non-stationary statistics that arise from ampli- 
tude modulation can be described by their spectral coher- 
ence properties (e.g. iBertolotti et al.l Il995h . The various 
types of wave coherence may be differentiated vi a auto- 
corre l ation and fluctuation spect r al est i mates (e.g. CordesI 
[19761: lEdwards fc Stappe7^l20ul 12001 Uenet fc Gi]|l2004h 
that are outside the scope of the current treatment. 
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4.3. Incoherent Sum 

When one or more incoherent sources of radiation are 
added together, the resulting covariance matrix of the to- 
tal ensemble mean Stokes parameters is the sum of the 
covariance matrices of the individual sources. For exam- 
ple, unpolarized noise adds a term to the diagonal of the 
covariance matrix, such that 

thereby reducing both the ellipticity of the distribution of 
the polarization vector and the multiple correlation be- 
tween total and polarized intensities. 

The incoherent addition of signals with different polar- 
izations, especially the superposition of orthogonally po- 
larized modes, also decreases the degree of polarization of 
the mean Stokes parameters. If th e modes are covariant 
(e.g. IMcK mnon fe Stinebri ng 1998), then the covariance 
matrix of the sum includes cross-covariance terms. For 
example, consider the incoherent sum of two sources de- 
scribed by Stokes parameters A and B. If the intensities of 
the two modes are correlated, then the resulting covariance 
matrix is given by 

C = C A + C S + S + S T , (32) 

where Ca and Cb are defined as in equation (|28[) and the 
cross-covariance matrix is 

= qsa^b AiBj , (33) 

where g is the intensity correlation coefficient. As shown 
in Appendix [Cj the incoherent superposition of covariant 
orthogonally polarized modes causes the variance of the 
total intensity to increase while that of the major polar- 
ization decreases. (The major polarization is defined by 
the eigen decomposition and is parallel to the major axis 
of the spheroidal distribution of S.) Furthermore, it is 
shown that the four-dimensional covariance matrix of the 
Stokes parameters can be used to derive the correlation 
coefficient as well as the intensities and degrees of polar- 
ization of superposed covariant orthogonal modes. 

4.4. Randomly Polarized Radiation 

iMcKinnonl (|2004h and lEdwards fe Stapperi (|2004l ) inde- 
pendently developed and applied novel eigenvalue analyses 
of the three-dimensional covariance matrix of the Stokes 
polarization vector S. Each noted an apparent excess dis- 
persion of S and concluded that it is due to the incoherent 
addition of randomly polarized radiation intrinsic to the 
pulsar signal. This hypothesis is based on the assumption 
that, apart from the proposed randomly polarized compo- 
nent, the noise in each of the Stokes parameters is purely 
instrumental, a premise that breaks down for sources as 
bright as the pu lsars studie d in th ese experiments. 

For example, IMcKinnonl |2004) analyzed observations 
of PSRs B2020+28 and B1929+10 that were recorded 
with the Arecibo 300 m antenna when the forward gain 
was 8 K/Jy and the system temperature was 40 K 
(Sti nebring et al.lll984h . Refer ring to the avera ge pulse 
profiles presented in Fig ure 2 of lMcK mno n1 (|2004h . the to- 
tal intensity of PSR B2020+28 peaks at 8 Jy where the 
modulation index is ~ 1.3. That is, the noise intrinsic to 
the pulsar signal exceeds the system equivalent flux den- 
si ty (SEFD; ~ 5 Jy) by as muc h as 100% (cf. Figure 6 
of [McKinnon fc Stinebrin"g1l2000t h Similarly, the self noise 



of PSR B1929+10 is as much as 50% of the SEFD. In 
both cases, source-intrinsic noise statistics cannot be ne- 
glected. To quantify the i mpact of self noise on single pulse 
polarimet ry the results of McKinnon (2004) are revisited. 
Note that lEdwards fc Stappersl (|2004h focused on the non- 
orthogonal modes of PSR B0329+54 and reported neither 
the instrumental sensitivity nor the statistics of the total 
intensity; therefore, this experiment i s not reviewed here. 

Figures 3 and 4 of IMcK mno i] (pool present two- 
dimensional projections of the ellipsoidal distributions of 
S at different pulse phases, along with the dimensions of 
the major axis a± and minor axes, a 2 and 0,3, at each phase. 
Panel a) in each of these plots indicates the off-pulse, in- 
strumental noise a n in each component of the Stokes po- 
larization vector; for PSR B2020+28, cr„ ~ 0.1 Jy and, for 
PSR B1929+10, a n ~ 0.04 Jy. Figures 2-4 are summa- 
rized in Table 1, which lists the pulse phase bin number; 
the modulation index the total intensity So; & n d the 
standard deviations of the total intensity o~q — /3So, the 
major polarization o\ = a\, and the minor polarizations 
v±_ = a 2 ~ a 3 . 

As shown in § 13.31 the standard deviation of the total 
intensity and the major axis of the spheroidal distribu- 
tion of the polarization vector should be equal; however, 
for every phase bin listed in Table 1, ao > o\. That is, 
there is no excess dispersion of the polarization vector and 
therefore no need for additional randomly polarized radi- 
ation; rather, the excess dispersion of the total intensity 
requires explanation. Noting that scalar amplitude modu- 
lation and wave coherence uniformly inflate the covariance 
matrix of the Stokes parameters, one possible explanation 
for Co > °~i is the incoherent superposition of covariant or- 
thogonally polarized modes (see Appendix ICj) . The coeffi- 
cient of correlation between the mode intensities g in equa- 
tions (|Cip and (|C2p is e quivalent to r\i in the analogou s 
equations (5) and (6) of iMcKinnon fc Stinebring (|1998[ ). 
The last term in all four of these equations indicates that 
positive correlation between the mode intensities increases 
the variance of the total intensity (<7q = Coo + ano - 
decreases that of the major polarization {a\ = C\\ + a^). 
(Anti-correlated mode intensities have the opposite effect.) 
The increased variance of the total intensity also explains 
why orthogonally polarized mode s typically coinci de with 
increased amplitude modulation ([McKin non 2001). 

Similarly, equation (|C3[1 shows that the self noise of 
partially polarized modes increases the dimensions of the 
minor axes (a| = C22 + c 2 and a\ = C33 + cr^). The 
source-intrinsic contribution diminishes to zero only for 
100% polarized modes. For all of the phase bins listed 
in Table 1, a± > a n , indicating that the modes are 
not completely polarized, as has been previously assumed 
(IMcKinnon fe Stinebringl 119981 : IMcKinnon fc Stinebringl 
2000). 

Equations (jCl|) to (jC4[) and the discussion in Ap- 
pendix O motivate the development of a new technique 
for producing mode-separated profiles. From the four- 
dimensional covariance matrix of the Stokes parameters, 
it is possible to derive the mode intensity correlation coef- 
ficient as well as the intensities and degrees of polarization 
of the modes. These values and the mean Stokes parame- 
ters, computed as a function of pulse phase, can be used to 
decompose the pulse profile into the separate contributions 
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Table 1 

Eigenvalue Analysis from IMcKinnonI (|2004T ) 



Bin 





So (Jy) cr 


o-i 




PSR B2020+28 


61 


0.45 


5.5 2.48 


1.14 


0.59 


74 


0.19 


2.8 0.54 


0.33 


0.30 


91 


0.56 


3.1 1.74 


0.91 


0.71 


PSR B1929+10 


90 


0.55 


0.68 0.38 


0.32 


0.09 


119 


0.68 


1.35 0.92 


0.63 


0.13 


125 


0.65 


1.00 0.65 


0.46 


0.10 



of the orthogonal modes. 

4.5. Giant Pulses 

The estimation of the degree of polarization of giant 
pulses is beset by fundamental limitations. On one hand, 
giant pulses remain unresolved at the highest time resolu- 
tions achieved to date. On the other hand, a sufficiently 
large number of independent samples must be averaged be- 
fore an accurate estimate of polarization is possible. Even 
if many time samples are averaged, those samples may be 
correlated due to scattering on inhomogeneities in the in- 
terstellar medium and/or the mean may be dominated by 
a single unresolved giant nanopulse that greatly reduces 
the effective de grees of freedom. 

For example. iPopov et al. (2006) presented a histogram 
of the degree of circular polarization of giant pulses from 
the Crab pulsar observed at 600 MHz. The histogram is 
interpreted as evidence that each giant pulse is the sum of 
~ 100 nanopulses, each with 100% circular polarization. 
However, referring to Figure [5J the width of the distribu- 
tion is also consistent with that of completely unpolarized 
radiation with < 2 degrees of freedom. Although the data 
were averaged over 256 samples to 32 fis resolution, the 
characteristic timescale of scattering in these observations 
was estimated to be 45 ±5 /is; therefore, the effective num- 
ber of degrees of freedom is of the order of unity. 

Small-number statistics also limit the conclusions that 
can be drawn from the correla tions between giant pulse 
spectra presented in Figure 12 o f lCordes et ail(|2004h . The 
asymptotic correlation coefficient at small lag was inter- 
preted as evidence that nanopulses are highly polarized. 
However, with only one estimate of the correlation coeffi- 
cient at a lag less than 0.1 seconds, the extrapolation to 
nanosecond resolution is questionable. Therefore, the con- 
straint on the degree of polarization at this timescale is of 
negligible significance. 

5. CONCLUSION 

A four-dimensional statistical description of polarization 
is presented that exploits the homomorphism between the 
Lorentz group and the transformation properties of the 



Stokes parameters. Within this framework, a generalized 
expression for the covariance matrix of the Stokes parame- 
ters is developed and applied to the analysis of single-pulse 
polarization fluctuations. The consideration of source- 
intrinsic noise renders randomly polarized radiation un- 
necessary, explains the coincidence between increased am- 
plitude modulation and orthogonal mode fluctuations, and 
indicates that orthogonally polarized modes are only par- 
tially polarized. Furthermore, the four-dimensional covari- 
ance matrix of the Stokes parameters enables estimation of 
the mode intensities, degrees of polarization, intensity cor- 
relation coefficient, and effective degrees of freedom of co- 
variant, orthogonal, partially polarized modes. Measured 
as a function of pulse phase, these parameters may be used 
to produce mode-separated profiles without any assump- 
tions about the intrinsic degree of mode polarization. 

The formalism is also used to derive the first and sec- 
ond moments of the degree of polarization as a function 
of the intrinsic degree of polarization and the number of 
degrees of freedom of the stochastic process. These are 
used to demonstrate that giant pulse polarimetry is fun- 
damentally limited by systematic bias due to insufficient 
statistical degrees of freedom. The discussions of ampli- 
tude modulation in § 14. H and wave coherence in ? !4.2l serve 
to illustrate the difficulties in defining an appropriate av- 
erage when the signal is unresolved. In this regard, the 
approach employed in this paper is complimentary to pre- 
vious analyses that make use of auto-correlation functions 
and fluctuation spectra to separately measure the statistics 
of the total and polarized intensities. Useful new results 
may be derived by extending the techniques employed in 
these works to include the cross-correlation terms that de- 
scribe the statistical dependences between the Stokes pa- 
rameters. 

I am grateful to J. P. Macquart for fruitful discussions 
and K. Stovall for technical assistance with the computer 
algebra system. The insightful comments of the referee led 
to significant improvements to the manuscript. M. Bailes, 
R. Bhat, J. Verbiest and M. Walker also provided helpful 
feedback on the text. 
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APPENDIX 

MARGINAL DISTRIBUTIONS OF THE MINOR STOKES PARAMETERS 

Equations ([19|) and (|20|) could have been also derived by first expressing equation (fl4|) in the natural basis. Here, it 
takes the simplified form, 

/(|eo|,|e 1 U) = ^-/o(|e |)/i(|ei|) ) (Al) 

where 

/ m (M) = ^exp(±^) (A2) 

are Rayleigh distributions. Note that in the natural basis, \eo\, |ei|, and tp are statistically independent and tp is uniformly 
distributed. The distributions of sq = |eo| + |ei| and si = |eo| — |ei| may then be obtained f r om th e convolution and 
cross-correlation, respectively, of /o(|eo|) and /i(|ei|). Similarly, following iFercher fc Steege"rl (|1981|) . equation (|A1|1 is 
used to compute the marginal distribution of s 2 = 2|eo||ei| cos-0. In the natural cylindrical coordinates defined in § 13.11 
the distribution of the radial dimension A = (s| + s §) -1 ^ 2 = 2|eo||ei| is found by integrating equation (|A1[) over ip, then 
transforming the integrand to yield the intermediate result 

/(A) = [°° foQeoDh f-£-)-L-d\e \ = -±-K (-£=) , (A3) 
Jo V 2 l e ol/ 2 l e o| A Ai VvAoAi/ 

where Kq is a modified Bessel function of the second kind. As the azimuthal dimension tp is uniformly distributed, the 
probability density of cos tp is 

/(cost/-) = ■ 

7ry 1 — cos z tp 

Combining with /(A) and performing another integral transformation yields 

' '' s 2 \ 1 , , 1 ( 2s 2 



/(s 2 )= / /(cost/')/ -dcosip = -exp[ — - (A4) 

Jo \ cos ip ) cos tp b \ b J 

for s 2 > 0, where S is the Lorentz interval. By symmetry, a similar result is found for s 2 < 0. 

DISTRIBUTION AND MOMENTS OF THE LOCAL MEAN DEGREE OF POLARIZATION 

Conversion of equation (|25j) to spherical coordinates and integration over all orientations of the polarization vector s 
yields the joint distribution of the local mean total and polarized intensities, 

, rrn 4n 2»-i|g|-2(»-2) / 2n3 S \ . { 2n\s\\S\ \ 

This joint density is defined on < \s\ < §o and is used to derive the distribution of the local mean degree of polarization, 
p = \s\/§q, as a function of the intrinsic degree of polarization P = \S\/Sq and the number of samples averaged n: 

r 2r(2n - 1) ( p2 ~ !)>(P 2 " i)"" 2 ([1 - P P^ 2n - t 1 + P P\^ 2n ) 

Hp)- I ^^■^°= r( ;_ 1 )r(„) -• ^ 

As P -> 0, 

= r f (2n ~ r )\ h^ 2 - 2n (^-l)pH P 2 - l)- 2 . (B3) 
1 (n — 1)1 (n) 

The distribution of the local mean degree of polarization is plotted in Figure [21 its first and second moments are 

(p) = J\f(p)dp=(l-P 2 rr(n+±) 3 F 2 (^2,n,n+l^,n+l;P 2 ^ (B4) 

and 

3(1 -P 2 ) n - /5 13 3 



(p 2 ) = / o P 2 / (^=^ TT ^ 3 ^^-,n,, + -;^n + ^;^j, (B5) 

where 3 F 2 is a regularized generalized hypergeometric function ([Wolfram Research! I2007D . These moments are used to 
calculate the variance of the local mean degree of polarization a 2 = (p 2 ) — (p) 2 . The theoretical values of (p) and a p are 
plotted as a function of n and P in Figures [3] and [4] 
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COVARIANT ORTHOGONAL PARTIALLY POLARIZED MODES 

Following the discussion in § 14.31 the covariance matrix of orthogonally polarized modes with correlated intensities is 
derived by starting with equation (f32|) . neglecting instrumental noise and considering only the self noise of the modes. If 
the modes A and B are orthogonally polarized, then A ■ B = — furthermore, if A is the dominant mode, then 

\A\ > \B\ and in the natural basis defined by A, the six non-zero elements of the covariance matrix are 

Coo - ?1PI| 2 + 4l|5|| 2 + 2gq A ^ B A Bo (CI) 

C 11 =<&\\A\\ 2 + <&\\B\\ 2 -2qs a s b \A\\B\ (C2) 

C 2 2 = C 33 =e A A 2 +e B B 2 (C3) 

Coi = Cio = 2^ 2 A A \A\ - 2&Bo\B\ + g^ B (B \A\ - A \B\). (C4) 

Including S — A + B and \S\ = \A\ — \B\, there are a total of seven unknowns and six unique constraints. However, if it 
is assumed that the modes have simi lar degrees of free dom (i.e. sa — Sb), the n the system can be solved numerically; e.g. 
using the Newton-Raphson method (jPress et al.lll992f) . In lMcK innonl (|2004l ). Cqi is not measured; therefore, no derived 
parameter estimates are currently presented. 
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